(submitted to Phys. Rev. B) 

Effect of a humid environment on the surface structure of RuO2(110) 
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Combining density-functional theory and thermodynamics we compute the phase diagram of 
surface structures of rutile RuO2(110) in equilibrium with water vapor in the complete range of 
experimentally accessible gas phase conditions. Through the formation of hydroxyl or water-like 
groups, already lowest concentrations of hydrogen in the gas phase are sufficient to stabilize an 
oxygen-rich polar oxide termination even at very low oxygen pressure. 
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I. INTRODUCTION 



Metal oxides are employed in a wide range of techno- 
logical applications, as for example catalysis (where they 
may serve as the support or the active material) and 
corrosion or wear protection in mechanical systems. A 
microscopic understanding of the function of these com- 
pounds requires knowledge of their often complex surface 
atomic structures, ^'^ which may be sensitively influenced 
by temperature and partial pressures in the surround- 
ing environment. This dependence could be particularly 
prominent in oxygen-containing environments, where ox- 
ide terminations of varying stoichiometry are easily per- 
ceived as a function of the abundance of oxygen in the 
surrounding. The situation becomes even more compli- 
cated by realizing that also hydrogen is almost always 
present in real- life systems, e.g., as a species in the reac- 
tant feed gas or as a bulk impurity. Then, it is in fact 
very plausible that surfaces are terminated by hydroxyl 
groups, as has recently been exemplified by Wang et al. 
for the Q;-Al2O3(0001) surface.^ 

Although it is quite obvious that such variations in the 
oxide surface structure and composition may affect the 
function of these materials in applications, systematic 
studies analyzing this issue are surprisingly sparse, even 
for well-defined oxide surfaces. ^'^ To this end we describe 
in the following how the method of first-principles atom- 
istic thermodynamics^^^'^ can be employed to construct 
a phase diagram of the most stable surface structures of 
an arbitrary oxide surface in equilibrium with its environ- 
ment in the complete range of experimentally accessible 
temperature and pressure conditions. As specific exam- 
ple we choose the RuO2(110) surface in contact with a 
humid gas phase formed of hydrogen and oxygen. Recent 
experiments have not only shown this surface to be an 
efficient CO oxidation catalyst but also a sensitive 
dependence of the measured turnover rate on the amount 
of hydrogen present in the system was reported^^. In 
these measurements, RuO2(110) showed a particularly 
high reactivity for the low temperature oxidation of CO 
by humid air, i.e. conditions where most technological 
catalysts cannot be applied due to anodic corrosion^^. 
This renders this oxide an interesting object of study de- 



spite its problematic technological applicability due to 
the poisonous nature of Ru. 

Similar to the results obtained by Wang et al. for 
Q-Al2O3(0001)^ we show in the following that also for 
RuO2(110) already lowest amounts of hydrogen in the 
gas phase lead to an overwhelming stability of hydro- 
genated surface phases. Through the formation of hy- 
droxyl or water-like groups an oxygen-rich polar oxide 
termination is then stabilized over the complete range of 
accessible oxygen pressures. Concomitantly, the strength 
of the 0-H bonds formed is so high, that even in ultra- 
high vacuum (UHV) annealing temperatures of the order 
of 600 K will be required to fully remove hydrogen from 
this transition-metal (TM) oxide surface. While this not 
only indicates the danger of unintentional hydrogen con- 
tamination in Surface Science experiments, our findings 
also demonstrate the sensitive dependence of oxide sur- 
face structure and composition on the surrounding gas 
phase. The latter should be systematically taken into 
account when addressing the functionality of these sur- 
faces in technical applications. 



II. THEORY 

In equilibrium with a given environment, characterized 
by temperature, T, and partial pressures, {pi}, the sta- 
ble surface structure minimizes the surface free energy, 
defined as 



Here, G{T,{pi],Ni) is the Gibbs free energy of the fi- 
nite crystal, Ni and ^i(T,pi) are number and chemical 
potential of the species of zth type, and 7(T, {pi\) is nor- 
malized to energy per unit area by dividing through the 
surface area A. In the present application to RuO2(110) 
in contact with a gas phase containing oxygen and hy- 
drogen, the surface free energy is thus a function of the 
three chemical potentials of Ru, O, and H. Assuming that 
the surface is also in equilibrium with the underlying ox- 
ide bulk constrains the chemical potentials of the two 
oxide components, as their stoichiometric sum has then 
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TABLE I: Experimental and theoretical binding energies, E\,, 
as well as vibration frequencies of isolated H2, O2 and H2O 
molecules. The theoretical values are obtained within our 

DFT-GGA approach using the same basis set that is also 
employed in the slab computations. 



H2 exp 

theory 
O2 exp 

theory 
H2O exp 

theory 



4.52 eV/mol.^^ 546 meV^'^ 
4.50 oV/mol. 538 meV 
5.16 eV/mol.^^ 193 meV^^ 
5.86 eV/mol. 180 meV 
9.60 eV/mol.^^ 454/466 mcV^'' 
9.81 eV/mol. 424/435 meV 



198 moV^'' 
189 meV 



FIG. 1: Vibrational contribution of surface hydrogen groups 
to the surface free energy of various hydrogenated phases of 
RuO2(110). See text for an explanation of the nomenclature 
used to describe the individual phases. 



to equal the Gibbs free energy of the bulk oxide. As a 
consequence the surface free energy can be discussed just 
in terms of a dependence on fio and ^h, which are both 
determined by the surrounding gas phase reservoirs and 
thus directly related to the experimentally controllable 
parameters, temperature and partial pressure®' 

We also note that not the absolute value of 
G{T, {pi}, Ni) enters the calculation of ^/{T, {pi}), but 
only the difference of the surface and the bulk Gibbs free 
energies. In ref. 9 we could show that for RuO2(110) in 
an oxygen environment the vibrational and entropic con- 
tributions to this difference, 7ru02(-^)' cancel to a large 
degree. This does, however, not comprise the contribu- 
tions from adsorbate functional groups. Since they have 
no counterpart in the bulk Gibbs free energy, their ab- 
solute vibrational energy and entropy enters into a con- 
tribution 7^''' (T) to the surface free energy of the cor- 
responding phase®. In most cases involving atomic or 
small molecular adsorbatcs this extra-contribution is for- 
tunately still negligible within our aspired accuracy of 
±10 — 15 meV/A^, as there is only a limited number of 
(low-energy) vibrational modes that gives rise to it. 

Due to their high-energy modes, hydrogen functional 
groups form a notable exception to this line of thought, 
yielding rather high 7'^'*^- (T) already for hydroxyl groups 
or water-like species at the surface. This is illustrated 
in Fig. 1 for various such phases discussed later in the 
text. The 7"*"^ (T) have been obtained from a detailed 
vibrational analysis of the various surface species, diago- 
nalizing the complete dynamic matrix while leaving the 
substrate fixedi^. We see that for T < TOOK 7^"^ (r) is 
of the order of 15 meV/A^ for hydroxyl groups and of 
the order of 30 meV/A^ for water-like adspccies. Inter- 
estingly, there is little variation in these values, when the 
functional groups adsorb in different surface sites, or even 
coadsorb closely to each other (see text below for the spe- 
cific nomenclature employed to describe the individual 



surface phases). Very similar 7^''^ (T) are also obtained 
if we use typical vibrational modes of these two surface 
functional groups when adsorbed on metallic surfaces^^. 
Although the individual modes change, these variations 
have apparently little effect on the composite quantity 
7^''^ (T). Furthermore, if we restrict our discussion to 
T < 700-800 K, also the temperature-dependence of this 
vibrational contribution to the surface free energy be- 
comes very small, cf. Fig. 1. Correspondingly, we may 
approximately take it into account by simply adding a 
constant 7^'^- of 15 meV/A^ per present hydroxyl-group 
and 30 me V/A2 per present water-like species in the com- 
putation of the 7(T, {pi}) of corresponding hydrogenated 
surface phases. 

Not aiming at an accuracy better than the stated 10-15 
meV/ A^, we may then replace the two free energies enter- 
ing 7(T, {pi}) by their total energies, and calculate them 
using density-functional theory within the full- potential 
linear augmented plane wave (FP-LAPW) approach^^ 
and the generalized gradient approximation (GGA) for 
the exchange-correlation functional^®. The high accu- 
racy of the basis set^° employed in these supercell calcu- 
lations has been detailed in a preceding publication that 
discussed the stability of RuO2(110) in a pure oxygen 
environment®. Thus, the numerical accuracy of differ- 
ences between j{T, {pi}) values of different surface struc- 
tures is better than ±5 meV/A^, while we stress that the 
aforedescribed approximate treatment of vibrational and 
entropic contributions docs not allow to distinguish be- 
tween phases which ^{T, {pi}) differ by less than about 
10-15 meV/A^. In comparison to our previous work we 
were forced to reduce the oxygen mufRn-tin radius to 1 
bohr, in order to treat the short 0-H bonds in the present 
extension to a gas phase containing both oxygen and hy- 
drogen. Consequently we had to increase the plane wave 
cutoff in the interstitial region from 17 Ry to 20.25 Ry 
for the wavefunctions and from 169 Ry to 400 Ry for the 
potential to maintain the same level of accuracy. Oth- 
erwise the same set-up is used for our computations as 
detailed before®. 

To further assess our ability to properly describe par- 
ticularly the energetics of O-H bonds within our DFT- 
GGA approach, we computed the binding energy, as well 
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FIG. 2: Top view of the RuO2(110) surface explaining the lo- 
cation of the two prominent adsorption sites (bridge and cus), 
which continue the bulk-stacking sequence (central panel). 
Additionally shown are perspective views of the two most sta- 
ble (1 X 1) terminations involving only O (left panel). The ge- 
ometry of hydrogenated surfaces is exemplified for two cases, 
comprising once a hydroxyl group at the bridge site (top right) 
and once adsorbed water at the cus site (bottom right). See 
text for explanation of the short hand notation used to de- 
scribe the various surface phases. Ru — light, large spheres, 
O = dark, medium spheres, H — white, small spheres. 



as vibration frequencies of free H2 and H2O molecules 
in the gas phase by diagonalizing their dynamic matrix. 
The results are compiled in Table I and compared to our 
results for molecular O2 . In contrast to the known strong 
overbinding in the latter case^^, the hydrogen-related re- 
sults show close agreement with reported experimental 
values. With a procedure described in detail in our pre- 
vious publication^, we can significantly reduce the large 
error introduced into the computed surface free energy 
due to the oxygen overbinding, giving us finally good 
confidence that the chosen set-up is adequate. Checking 
especially on the approximate nature of the exchange- 
correlation functional we additionally repeated the most 
relevant computations entering into the construction of 
the surface phase diagram within the local-density ap- 
proximation (LDA)^** and will comment on these findings 
in the discussion below. 



III. SURFACE STRUCTURE OF RUOzCllO) 

RuO2(110) has recently attracted considerable inter- 
est as playing a key role in the catalytic CO oxidation 
at a Ru(OOOl) surface that is exposed to high oxygen 
pressure. The atomic arrangement of this ru- 
tile (110) surface, also known from TiO2(110), exhibits 
two distinct adsorption sites within the computed (6.43 
A X 3.12 A) rectangular surface unit-cell. Figure 2 shows 
these two prominent sites over the oxygen-poor termina- 
tion of RuO2(110): The so-called coordinatively unsatu- 
rated (cus) site (atop of fivefold coordinated Ru atoms, 
Ru'^"''), and the bridge (br) site between two fourfold co- 
ordinated Ru atoms (Ru'^''). The bulk stacking sequence 
would be continued by oxygen atoms first occupying the 
bridge (O'^'') and then the cus (O'^"'') sites, thus leading 
to the two other (1 x 1) terminations of a rutile (110) 



surface shown in perspective view in the left panels of 
Fig. 2, cf. ref. 9 for more details. 



To evaluate the surface phase diagram we computed 
what we believe are all possible (1 x 1) (co-)adsorption 
geometries of O and H at the cus and bridge sites of 
RuO2(110). Additionally, we also considered any pos- 
sible combination of adsorbed molecular H2, OH, and 
H2O groups at these two adsorption sites. In the fol- 
lowing we will use a short hand notation to characterize 
this manifold of studied geometries, indicating first the 
occupancy of the bridge site and then of the cus site, e.g. 
{OH)^" /— for an OH-group at the bridge site, the cus 
site being empty, cf. Fig. 2. Although adsorption of H 
or O is in principle also conceivable at other sites on the 
surface, test calculations always resulted in significantly 
lower binding energies compared to adsorption at the cus 
or bridge sites. Correspondingly, we may restrict our dis- 
cussion to these two prominent sites in the following. In 
view of the more than 20 studied geometries we would 
like to emphasize at this point that a complete structural 
relaxation - including symmetry breaking at the surface 
- was essential to obtain the correct energetics, as we 
find particularly the tilting of OH and water-like groups 
to decisively influence the relative stabilities of the tested 
geometries. 



Our self-consistent DFT total energies implicitly con- 
tain all lateral interactions within the considered (1x1) 
surface unit-cell. Significant lateral interactions beyond 
the first-neighbor cus-br interaction, cf. Fig. 2, could e.g. 
lead to ordered adsorption phases that are more dilute. 
For pure oxygen phases we therefore performed a sys- 
tematic study to extract all lateral interactions between 
adsorbed oxygen atoms employing large surface unit-cells 
up to (4 X 1) and (2 x 2) periodicity^^. We found only 
very small lateral interactions, which we understand be- 
ing a consequence of the rather open surface structure 
(compared to the more often studied close-packed metal 
surfaces). We don't expect this picture to change sig- 
nificantly in the case of H and O coadsorption, as con- 
firmed by a series of test calculations employing (2x1) 
surface unit-cells. The remaining small lateral interac- 
tions or disorder effects will only affect the relative sta- 
bilities close to the boundaries between different stable 
phases. In the present study we aim to construct the 
large scale phase diagram in the whole range of exper- 
imentally accessible gas phase parameters. Small vari- 
ations in the vicinity of the obtained phase boundaries 
do not change the physical conclusions we deduce on the 
basis of this phase diagram and are thus not of present 
concern. Hence, we will restrict the discussion here to the 
(still large) number of (1 x 1) geometries and defer a sys- 
tematic study of larger unit-cell phases to a forthcoming 
publication. 



4 



IV. SURFACE PHASE DIAGRAM OF 
HYDROGENATED RUO2(110) 

When now discussing the dependence of the computed 
surface free energies on the chemical potentials of the gas 
phase species, we note that experimentally (and assum- 
ing that thermodynamic equilibrium applies) /io cannot 
be varied without bounds. If it gets too low, all lattice 
oxygen atoms would preferably go into the gas phase, 
i.e the oxide would decompose. On the other hand, if 
Ho gets too high, oxygen will start to condense on the 
surface. Theoretical estimates for both extremes can be 
derived^ and we will correspondingly restrict the presen- 
tation of the surface free energies to the limited range 
between this O-poor and 0-rich boundary, respectively. 

Extending this to a multi-component environment we 
note that depending on the gas phase reaction barrier (s) 
an equilibrium composition may or may not be present in 
the gas phase. If the barrier is rather high, as in the case 
of the CO oxidation reaction, the first-principles atom- 
istic thermodynamics concept can e.g. be employed to 
address the situation of a constrained equilibrium where 
each gas phase species is separately in equilibrium with 
the surface, but not with the respective other gas compo- 
nent (s). Depending on the reaction barriers for the ad- 
sorbed species, surface kinetic effects may then become 
important for certain temperature and pressure condi- 
tions, so that the actual surface structure deviates from 
the prediction of the thermodynamic approach^". 

This problem does not arise, if one concentrates on 
the situation where the environment is itself fully equili- 
brated. In the present context of a gas phase containing 
oxygen and hydrogen this corresponds to allowing for the 
possibility of water formation in the oxyhydrogen reac- 
tion. In fact, over the whole finite stability range possible 
over RuO2(110), /io is already so high that water is not 
only a possibility, but the dominant gas phase compo- 
nent, if such a full equilibrium is reached^. Modelling this 
particular situation, the chemical potential of hydrogen 
atoms at the RuO2(110) surface is therefore determined 
by a H2O reservoir in the following. To assure that we 
restrict the discussion to a humid environment, i.e. stay 
within the gaseous state of water for any temperature or 
pressure, the maximum /XH2O considered will be at -0.91 
eV with respect to the total energy of a H2O molecule, 
which corresponds to the chemical potential of water at 
the (experimental) critical point. 

Before addressing the equilibrium with the two com- 
ponent gas phase, it is instructive to first discuss the 
oxide surface structure in a pure O2 atmosphere^. Then, 
the surface free energy depends only on the chemi- 
cal potential of oxygen, which in turn represents the 
two-dimensional dependence on temperature and oxygen 
pressure. Figure 3 shows the limited stability range be- 
tween the O-poor and 0-rich boundary. To make the de- 
pendence on /io a bit more lucid, we additionally show 
a pressure scale in the top x-axis for fixed T = 600 K, 
which is a typical annealing temperature employed in 
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FIG. 3: Surface free energies of the three (1 x l)-RuO2(110) 
terminations involving only oxygen (solid lines). The corre- 
sponding 7(T, {pi}) of hydrogen containing phases are drawn 
with dashed lines, if no O is present at the cus sites, and with 
dotted Hnes, if O is present (see Fig. 2 for nomenclature and 
site explanation). The dotted vertical lines indicate the al- 
lowed range of oxygen chemical potential, while the chemical 
potential of H is in this plot chosen for very H-rich condi- 
tions corresponding to the experimental critical point of wa- 
ter (/tHaO = — O.OleV'*). In the top x-axis, the dependence 
on /io (T, PO2 ) has been cast into a pressure scale at a fixed 
temperature of T = 600K. 



this system^^'^^. The solid lines indicate the surface free 
energies of the three possible (1x1) terminations of a 
rutilc (110) surface. While the mixed 0-Ru termination, 
where both cus and bridge sites are unoccupied (-/-), is 
rather high up in energy, the two lines corresponding to 
the other two terminations displayed in Fig. 3 intersect in 
about the middle of the possible range of oxygen chemical 
potential. At lower /io, corresponding to ultrahigh vac- 
uum (UHV) conditions, the traditionally discussed stoi- 
chiometric termination is stable, where oxygen occupies 
only the bridge sites at the surface (O^''/— ). However, 
towards higher oxygen pressures we find a second, 0-rich 
termination to become more stable, in which O"^"® atoms 
arc now also present (O'^'^/O'^"''). 

When next addressing the equilibrium with a gas 
phase containing additionally water vapor, the one- 
dimensional dependence of the surface free energy on 
/io gets extended to the two-dimensional dependence on 
(McMHso)- The maximum effect of the additional gas 
component on the surface structure can be seen from the 
surface free energies of hydrogenated geometries in the 
upper H-rich linut of /iHaO described above"'^^. For this 
particular chemical potential (/iHaO = -0.91 eV) the lines 
corresponding to the most stable structures studied are 
included in Fig. 3: Dotted lines represent geometries 
including O'^"*' atoms, while dashed lines depict phases, 
where no oxygen is present at the cus sites. Interestingly, 
we find several cases, where different surface geometries 
of identical stoichiometry, which correspondingly exhibit 




FIG. 4: Surface free energies of the four stable phases (or 
degenerate phase pairs) discussed in the text, as a function 
of (/io,MH2o)- The corresponding phases are drawn as thick 
Unes in Fig. 3 for the H-rich Umit. 



parallel j{T, {pi}) linos in Fig. 3. turn ont energeti- 
cally virtually degenerate, i.e. their surface free energies 
within DFT-GGA differ by less than ±5 meV/A^ and are 
thus not distinguishable within the approximations of our 
approach described above. In such cases we only include 
one line in the graph and label it with both geometries. 

Particularly relevant for the future discussion will be 
the two thicker drawn dotted lines, each representing a 
degenerate pair of surface phases. In the first case, two 
H atoms are present per (1 x 1) surface unit-cell, forming 
either a chemisorbed water-like species at the cus sites 
(0^'"/(H20)'^"'^) or hydroxyl groups at both bridge and 
cus sites ((0H)'^7(0H)'="''). In the second case, one H 
atom per unit-cell forms a hydroxyl group either at the 
bridge or at the cus site ((0H)''70'="" and 0''7(0H)™" 
respectively), cf. Fig. 2. We note that particularly the 
virtual degeneracy of these last two phases is somewhat 
surprising, as they correspond to quite different geome- 
tries, namely to bridge-bonded and atop-bonded OH- 
groups at different surface sites. As apparent from Fig. 3, 
the surface free energy in this H-rich limit is significantly 
lowered in both cases (thick dotted lines) in comparison 
to the two stable pure O terminations (thick solid lines) 
discussed before. 

Allowing variations of /iHaO away from the H-rich limit, 
the two-dimensional graph in Fig. 3 turns into the; thrc;e- 
dimensional one shown in Fig. 4, where we now only 
include the surface free energy planes of the four most 
stable phases (or degenerate phase pairs), that will play 
a role in the later discussion. Note that at the maximum 
MH2O = -0.91 eV along the front a;-axis in Fig. 4 we 
simply recover the two-dimensional dependence of Fig. 
3. While the j{T, {pi}) of the two pure O terminations, 
0^"/~ and O^'^/O™^, do, of course, not depend on fin^O, 
the hydrogenated geometries become increasingly less fa- 
vorable the lower the water chemical potential becomes, 
i.e. the less hydrogen is present in the gas phase. For the 
lowest /iHaO shown in Fig. 4, all H involving phases have 
then already become less stable than the pure O termi- 
nations for any /io, i.e. the H content in the gas phase is 
so low, that no hydrogen can be stabilized at the surface 
anymore. 

Although instructive, the three-dimensional presenta- 
tion does not allow an easy access to the most relevant 
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FIG. 5: Surface phase diagram of stable structures of 
RuO2(110) in equilibrium with a humid environment, i.e. as 
a function of (/io,/iH2o) in the gas phase. The additional 
axes show the corresponding pressure scales at T = 300 K 
and 600 K. 



information contained in Fig. 4, namely which phase has 
the lowest surface free energy for given gas-phase condi- 
tions, i.e. given (/io, fJ-a^o)- Therefore we display in Fig. 
5 this surface phase diagram of stable structures, which 
corresponds to taking a look at Fig. 4 from below. We 
find all four phases (or degenerate phase pairs) discussed 
up to now to be most stable in different (/^Qj /-''H20)- 
regions: At the lowest fin^o shown so little hydrogen is 
present in the environment, that just the surface struc- 
ture dependence of RuO2(110) in a pure O2 atmosphere 
is recovered, i.e. the stoichiometric O'"'/— phase is stable 
at low /io, while towards higher oxygen pressures, O™'' 
is additionally stabilized leading to the o'^''/0'^"'' phase. 
Upon increasing the H content in the gas phase hydro- 
gen first forms one hydroxyl group per surface unit-cell 
((0H)''70™" or 0''7(0Hj™''), whereas towards H-rich 
and 0-poor conditions in the upper left part of Fig. 5 
the maximum coverage of two H atoms per surface unit- 
cell is reached. This leads either to two hydroxyl groups 
at both bridge and cus sites ( (0H)^7(0H)"''' ) or to a 
chemisorbed water-like species, bound with leV to the 
cus sites as shown in the right bottom panel of Fig. 2 
(0''7(H20)^"«). 

In order to check how much our results may be affected 
by the exchange-correlation functional, we recomputed 
the stable surface phases using the LDA, i.e. consider- 
ing the LDA lattice constant for the substrate (6.45 A 
X 3.16 A) and including a full geometry reoptimization. 
Although the absolute surface free energies are found to 
be on average about « 15meV/A^ higher compared to 
the DFT-GGA results, the same energetic sequence as 
displayed in Fig. 3 was obtained. In particular, the 
two pairs of degenerate surface structures discussed be- 
fore give again surface free energies differing by less than 
15 meV/A^ from each other and thus exhibit stabilities 
that cannot be distinguished within the current approx- 
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imations. Correspondingly, the phase diagram obtained 
within the LDA exhibits the same four stable regions and 
the same qualitative structure as the one shown in Fig. 
5. Quantitatively, the phase boundaries of the hydro- 
genated phases are shifted by about 0.5 eV towards lower 
/iHaOj indicating that hydrogen bonds are stronger within 
the LDA. In line with our preceding study^ we there- 
fore conclude that the uncertainty in the DFT exchange- 
correlation functional may affect transition temperatures 
or pressures between stable phases deduced on the basis 
of the constructed phase diagrams. Yet the transitions 
per se, i.e. the overall structure of the phase diagram, 
are not changed. 

To translate the dependence on the chemical poten- 
tials into temperature and pressure, the pressure scales 
at room temperature and T = 600 K are also included in 
Fig. 5. From those we can derive e.g. what pressures 
would be required in the gas phase in order to stabilize 
hydrogen at RuO2(110). The transition from the pure 
O terminations to the H involving phases at somewhat 
academic hydrogen pressures of the order of 10~^° atm 
for T = 300 K simply indicates that the oxide surface will 
always be hydrogenated at this low temperature in any 
realistic environment, or in other words that the other 
areas of the phase diagram are experimentally not acces- 
sible at T=300K. From our results we can further con- 
clude, that even under UHV conditions, temperatures of 
the order of 500-600 K (GGA, or 600-700 K within the 
LDA) would be required to remove hydrogen completely 
from the surface, while at atmospheric pressure this rises 
to about 1000 K, i.e. to temperatures already close to 
the decomposition of the whole oxide. It is interesting to 
notice that in a number of recent experimental studies on 
RuO2(110) an annealing temperature of only T = 600 K 
was employed^^'^^'^^. According to our results this might 
not have been sufficient to completely remove residual H 
contaminants from the surface. Similarly, our results also 
imply that hydrogen has to be taken into account as a 
likely surface species under all realistic gas phase condi- 
tions up to rather elevated temperatures, e.g. also in all 
catalytically-relevant environments. 



V. COMPARISON WITH HYDROGENATED a- 
AI2O3 

The phase diagram in Fig. 5 reveals that the most no- 
table effect of hydrogenation is to stabilize O"^"^ atoms 
over the complete range of possible oxygen pressure: 
While in the pure O2 atmosphere the phase including 
O'^"'' becomes only more stable towards higher oxygen 
chemical potentials, this species is present in all four sta- 
ble hydrogenated phases for any /xq at sufficiently high 
MH2OJ cf. Fig. 5. In other words, already lowest concen- 
trations of hydrogen in the gas phase are sufficient to sta- 
bilize a non-stoichiometric surface terminated by either 
hydroxyl groups (one H at O) or water-like species (two 
H at O). This finding is analogous to the results on a- 



Al2O3(0001) obtained by Wang et al.^ , who reported the 
purely hydrogen-induced stability of the 0-rich 0-plane 
termination of this corundum structure. Yet, while in 
the latter case the polar termination is never favorable 
without the presence of hydrogen, a humid environment 
only increases the stability range of the 0-rich surface (in 
terms of hydroxyl groups) in the case of RuO2(II0). 

The hydrogenation goes in both cases hand in hand 
with a drastic reduction of the work function. For 
RuO2(II0) adding the C^"*^ atoms first significantly 
increases the work function by -1-1.5 eV, reflecting 
the increased dipole moment of the 0-rich o'"'/0'^"*' 
termination^. On the contrary, adding hydrogen into 
the surface geometry markedly reduces the work func- 
tion again by -2.0 eV or -3.5 eV for the phases involving 
one or two hydroxyl groups, respectively. On a-Al203, 
which was found to stabilize up to three hydroxyl groups 
per surface unit-cell, even larger work function reductions 
by more than 6eV were obtained^. 

Unfortunately it is not clear, whether the latter 3H 
per cell coverage corresponds already to full saturation, 
as the possibility of water-like species was not considered 
in the study by Wang et al..^ Nevertheless, already for 
the fully hydroxylated a-Al2O3(0001) surface negative 
surface free energies were obtained in the H-rich limit, 
indicating the strength of the H-OAl bond and consis- 
tent with the negative heat of the q;-A1203 -I- 3H2O — > 
2A1(0H)3 reaction'^. On the other hand, the weaker H- 
ORu bond never leads to negative 7(T, {pi}) even at the 
most H-rich conditions, cf. Fig. 3. This lower bond 
strength is also reflected by the result that the fully hy- 
droxylated RuO2(110) surface still exhibits a compres- 
sion of the topmost layer distance, although this compres- 
sion is noticeably reduced compared to the hydrogen free 
surface. In contrast it had been found that the topmost 
layer is already strongly relaxed outwards at the fully hy- 
droxylated corundum Al2O3(0001) surface^. Finally we 
note in passing that both studies interestingly obtain a 
rather low site specificity for the formation of OH-groups 
on the 0-rich surface, as in our case expressed by the 
degenerate surface free energies of the (OII)'^''/0'^"*' and 
0'""/(0H)'^"** phases. Recalling the large difference in 
binding energy between O'"' and O'^"'*, 2.5 eV vs. 1.2 eV 
with respect to molecular oxygen^'^", this result is rather 
unexpected. 



VI. SUMMARY 

In conclusion we presented the ab-initio based surface 
phase diagram of stable structures of RuO2(110) in equi- 
librium with a humid environment formed of oxygen and 
water vapor. We find the oxide surface to be completely 
hydrogenated at room temperature for any realistic gas 
phase conditions. Depending on the partial pressures in 
the gas phase H is either present in form of hydroxyl 
groups at bridge and/or cus sites, or in form of a water- 
like species at the cus sites of RuO2(110). A complete 
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structural relaxation allowing also for a tilting of the var- 
ious surface groups was found necessary in order to arrive 
at the correct relative stabilities of the various phases. 

The most notable effect of hydrogen on the oxide sur- 
face structure is the stabilization of a polar, 0-enriched 
surface structure through a considerable reduction of the 
finite dipole moment introduced upon addition of the cor- 
responding terminal oxygen atoms. This finding is analo- 
gous to the one of a previous study on the effect of hydro- 
genation on the stability of the a-Al2O3(0001) corundum 
surface. The lower H-ORu bond strength does however 
not lead to an instability of the bulk oxide structure, as 
reflected in the reported negative surface free energies of 
fully hydroxylated q;-A12O3(0001) surfaces under H-rich 
conditions. Still, the 0-H bonds formed on RuO2(110) 
are strong enough that even in UHV annealing temper- 
atures of the order of ~ 600 K will be required to com- 
pletely remove residual H contaminants from the surface. 



For atmospheric gas phase pressure this stability extends 
roughly up to 1000 K, indicating the necessity to include 
H as a likely surface species when modeling the function 
of Ru02 oxide surfaces in technological applications like 
e.g. catalysis. 
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